DESY 98-051 

LUNFD6/(NFFL-7156) 1998 



ISSN 0418-9833 



Resolved photon processes in DIS and small x dynamics. 



H. Jung, L. Jonsson, H. Kiister 
Physics Department, Lund University, P.O. Box 118, 221 00 Lund, Sweden 

Abstract 

It has been found that recent results on forward jet production from deep inelastic scat- 
tering can neither be reproduced by models which are based on leading order a s QCD 
matrix elements and parton showers nor by next-to-leading order calculations. The mea- 
surement of forward jet cross sections has been suggested as a promising probe of new 
small x dynamics and the question is whether these data provide an indication of this. 
The same question arises for other experimental data in deep inelastic scattering at small 
x which can not be described by conventional models for deep inelastic scattering. In this 
paper the influence of resolved photon processes has been investigated and it has been 
studied to what extent such processes are able to reproduce the data. It is shown that 
two DGLAP evolution chains from the hard scattering process towards the proton and 
the photon, respectively, are sufficient to describe effects, observed in the HERA data, 
which have been attributed to BFKL dynamics. 

1 Introduction 

The cross section of forward jet production in deep inelastic scattering (DIS) has been advo- 
cated as a particularly sensitive measure of small x parton dynamics [ [T], |2| . Analytic calcula- 
tions based on the BFKL equation in the leading logarithmic approximation (LLA) are in fair 
agreement with data. However, recent calculations of the BFKL kernel in the next-to-leading 
logarithmic approximation (NLLA) [ [| have given surprisingly large corrections, and it remains 
to be shown whether the data can still be reasonably described. 

Monte Carlo generators based on direct, point-like photon interactions (DIR model), cal- 
culated from leading order (order a s ) QCD matrix elements, and leading log parton showers 
based on the DGLAP evolution do not take any new parton dynamics in the small x region 
into account and are therefore not expected to fit the experimental data. Recent results from 
the HI [ U and ZEUS [ || experiments on forward jet production exhibit significant deviations 
from the predictions of such models. Also next-to-leading order calculations (NLO i.e. order 
a,) predict too small a cross section compared to data. 

Similarly the DIR model is unable to reproduce the fractional di-jet rate, R2, as a function 
of x and Q 2 [ §], investigated in a kinematic domain where Vr,jet ~ Q 2 ■ m these events the two 



high pt jets originate from the hard scattering process. Discrepancies up to a factor of three 
are observed. Inclusive jet cross sections in the transition region between photo-production and 
DIS are measured to be significantly larger than predicted by models assuming DGLAP based 
direct interactions [|7|]. 

Measurements of the transverse energy flow [ ||, [J do not provide a clear distinction between 
various models, whereas it has been demonstrated that the data on transverse momentum 
spectra of charged particles [ |10J can not be described by conventional leading order DGLAP 
equations at small x and Q 2 , a region where new parton dynamics are expected to become 
important. 

On the other hand all these data samples are rather well described by the color dipole model 
(CDM) as implemented in the Monte Carlo program ARIADNE In the CDM the gluon 

emission is described as radiation from color dipoles stretched between quark-antiquark pairs. 
These color dipoles radiate independently and therefore the gluons are not ordered in transverse 
momentum (fey). This dis-ordering in kj- is also a prominent feature of the BFKL dynamics 
where the parton emission follows a "random walk" in kx- The DGLAP evolution used in the 
DIR model gives strong ordering in kr- The question is thus whether all these data provide an 
indication of new small x parton dynamics. 

In this work we have studied the possible contribution of virtual resolved photons to the cross 
section of the DIR model as a possible alternative explanation for the differences between data 
and the Monte Carlo predictions. All the studies have been performed with the RAPGAP 2.06 [ 
0, |H| Monte Carlo event generator. It has been checked that this generator agrees well with 
analytical LO calculations and with the LEPTO 6.5 [ [TJ[] Monte Carlo in the case of direct 
interactions. 



2 Resolved Photons in DIS 

In electron-proton scattering the internal structure of the proton as well as of the exchanged 
photon can be resolved provided the scale of the hard subprocess is larger than the inverse 
radius of the proton, 1/R 2 ~ ^ 2 qcdi an d the photon, 1/R 2 ~ Q 2 , respectively. Resolved 
photon processes play an important role in photo-production of high jets, where Q 2 ~ 0, 
but they can also give considerable contributions to DIS processes [ |7|, |15[ if the scale /i 2 of the 
hard subprocess is larger than Q 2 , the inverse size of the photon. 

In the following we give a brief description of the model for resolved virtual photons used in 
the Monte Carlo generator RAPGAP. Given the fractional momentum transfer of the incoming 
electron to the exchanged photon, the Equivalent Photon Approximation provides the flux of 
virtual transversely polarized photons [ 0, |1^, and references therein]. The contribution from 
longitudinally polarized photons has been neglected. The partonic structure of the virtual 
photon is defined by parameterizations of the parton densities, x 7 / 7 (x 7 , fi 2 , Q 2 ), which depend 
on the two scales fi 2 and Q 2 [ [17|, fji| . The following hard subprocesses are considered (RES 



model): gg — ► qq, gg — > gg, qg — > qg, qq — > gg, qq — > qq, qq — > qq. Parton showers on both 
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Figure 1: Deep inelastic scattering with a resolved virtual photon and the q^g p — > qg partonic 
subprocess. 

the proton and the photon side are included. The generic diagram for the process q^g p — > qg 
including parton showers is shown in Fig. [I]. 

Since the photon structure function depends on the scale, fi 2 , of the hard scattering process, 
the cross section of resolved photon processes will consequently also depend on the choice of 
this scale. The parton density of the photon is evolved from a starting scale Ql to the scale /i 2 , 
the virtuality at the hard subprocess, giving a resummation to all orders. 

It should be noted that a NLO calculation assuming point-like virtual photons contains a 
significant part of what is attributed to the resolved structure of the virtual photon in the RES 
model [ [HJ. This is due to the fact that in order a 2 the virtual photon can split into a qq pair 
with one of the two quarks interacting with a parton from the proton, giving rise to two high 
Pt jets. 

2.1 Choice of Scale 

In leading order a s processes the renormalization scale and factorization scale [ip are not 
well defined which allows a number of reasonable choices. There are essentially two competing 
effects: a large scale suppresses a s (/z 2 ) but gives, on the other hand, an increased parton density, 
xf(x,fi 2 ), for a fixed small x value. The net effect depends on the details of the interaction 
and on the parton density parameterization. 
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However, in resolved virtual photon processes the choice of the scale /i 2 , at which the photon 
is probed, is severely restricted [ |20| . If we proceed from the electron (see Fig. [I]) towards the 
hard subprocess, the first branching occurs already at the electron vertex e — > e ,r y*, where the 
scale in DIS is chosen to be Q 2 = —t (t being the Mandelstam variable of a 2 — > 2 process), 
which is related to the transverse momentum of the scattered electron by: Q 2 = p\. e ,j{\ — y). 
Here y is defined by y = jp, with / and P being the four vectors of the incoming electron and 
proton and q is the four vector of the (virtual) photon with q 2 = —Q 2 . This indicates that the 
scale must be larger than p^. e ,, for a typical value of y ~ 0.5, this gives Q 2 = 2 • p\. e . 

Within a consistent picture of the evolution from the electron towards the hard scattering 
process, the definition of the scale should be kept the same for all branchings. Thus the proper 
scale for a hard scattering process in the t - channel (as shown in Fig. [I]) would be t, since 
Q 2 = —t. Here t refers to the hard scattering process. However not only t - channel diagrams 
contribute, but also u and s - channel diagrams as well as interference terms. The calculation 
of the relative contributions to the cross section in principally limited due to the interference 
terms. We are therefore forced to make a best guess of a reasonable scale. We have chosen 2 
different scales: fi 2 = 4p^ and fi 2 = Q 2 + p\- The first one is closely related to the t - channel 
diagrams, where we have included a factor of 4 to account for the fact that \t\ > p\. The 
second one is chosen for a smooth transition from usual DIS to DIS including resolved photons 
and to photo-production. This choice of scale has also been used in NLO calculations including 



resolved photons in deep inelastic scattering [ IS . 



2.2 Parton Distribution Functions 

Two parameterizations of the parton distribution in the proton, GRV 94 HO (DIS) and CTEQ4M 
have been considered, which both give good agreement with the proton structure function data 
[ ^T], p2] . It was found that the produced results were identical within the percent level, when 
keeping A QC d fixed. 

Also for the virtual photon two different parameterizations of the parton density have been 
investigated. The photon can interact via its partons either in a bound vector meson state or 
as decoupled partons if the pt of the partons is high enough. Only the latter is relevant in the 
Q 2 range considered in this paper. The splitting 7 — > qq is called the anomalous component of 
the photon. The SaS [ [17] parameterization offers a choice of Qq values at which the anomalous 



part becomes effective. We have studied these choices resulting in different magnitudes of the 
parton densities, and consequently of the cross sections. For the SaS parameterization we have 
used Ql as given by eq.(12) of ref. [ |17| (IP2 = 2), which gave the largest resolved photon 
contribution. 

The GRV LO parton density description of the real photon together with the virtual photon 
suppression factor of Drees and Godbole (DG) [ [0J gives similar results. Drees and Godbole 
proposed a x 7 - independent suppression factor: 
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with an adjustable parameter P 2 , which is set to P 2 = 0.5 GeV 2 typical for a hadronic scale. 
Comparisons between the SaS and the DG parameterizations are given for some of the observ- 
ables. 

The GRS [ structure function for the virtual photon has not been used, since it is 
restricted to Q 2 < 10 GeV 2 which is not useful for the studies presented here, 

If not stated otherwise, in the following we have used the scale /i 2 = Q 2 + p%, the CTEQ4M 
parameterization of the parton density in the proton and the DG description of the photon. 



3 Forward Jets 

HERA has extended the available x region down to values below 10~ 4 where new parton dy- 
namics might show up. Based on calculations in the LLA of the BFKL kernel, the cross section 
for DIS events at low x and large Q 2 with a high p\ jet in the proton direction (a forward 
jet) [ [l], |2j is expected to rise more rapidly with decreasing x than expected from DGLAP 
based calculations. New preliminary results from the HI [ || and ZEUS [ [5] experiments have 
recently been presented. The data can be described neither by conventional DIR Monte Carlo 
models nor by a NLO calculation, while comparisons to analytic calculations of the LLA BFKL 
mechanism has proven reasonable agreement. 

It should be kept in mind that both the NLO calculations and the BFKL based calculations 
are performed on the parton level whereas the data are on the level of hadrons. Nevertheless 
these comparisons make sense, considering the present measuring errors and the fact that 
simulations have shown that the hadronization effects are at most 20%. 

In this study, Monte Carlo events from the RAPGAP generator have been used to investigate 
whether the experimental data can be equally well reproduced by the inclusion of resolved 
photon processes. The cuts applied in the laboratory frame are equivalent to those in the 
preliminary HI analysis [ffl. 



Kinematic cuts: 



Jet selection: 
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where 8 = 0° corresponds to the proton beam direction. 



Jets were reconstructed in the rj-(f> space by applying a cone algorithm with a cone radius 
of R = ^/rj 1 + 2 = 1. The jet cuts have been defined to ensure that the momentum fraction of 
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the jet {xj et = Ejet/Ep, where E p is the proton energy) is large compared to x so as to maximize 
the phase space for a evolution in x. At large Xj et the parton distributions of the proton are well 
measured and the uncertainty from non-perturbative effects is avoided. In order to suppress 
the phase space for DGLAP evolution, the p\ of the jets was further required to be of the same 
order of magnitude as Q 2 (0.5 < p\IQ 2 < 2). Thus DIS events with small x and an energetic 
high pt jet in the forward direction are selected. 
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Figure 2: The forward jet cross section as a function of x using the cuts specified in the text. 
The data points are preliminary HI data [ffl. In a. are shown the RAPGAP predictions for 
the sum of direct and resolved processes using again the scales /i 2 = Q 2 + p\ (solid line) and 
fj 2 = 4p2, (dotted line). Also included are the predictions for direct DIS processes with the 
scales fi 2 = Q 2 + p\ (dashed line) and /i 2 = Ap\ (dash-dotted line), respectively. The parton 
distribution of the virtual photon was parameterized according to SaS 2 DIS. In b. the sum of 
direct and resolved processes (solid line) and the direct (dashed line) are presented using the 
DG + GRV description of the parton density in the photon and with the scale /x 2 = Q 2 + p\. 
In all cases the parton density of the proton was described by CTEQ4M. 

In Fig. the preliminary data of the HI collaboration [ |j] are compared to the RAPGAP 
Monte Carlo predictions for both, direct process alone (labeled DIR) and for the sum of direct 
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and resolved processes (labeled DIR+RES). The comparison is performed for the two scales, 
ji 2 = Q 2 +Pt an d = ^Pt an d for the two parameterizations of the photon parton density (SaS, 
Fig. ^|a, and DG, Fig. Qb). It is clearly seen that the direct contribution alone is not sufficient 
to reproduce the data whereas the DIR+RES model gives good agreement with data. The 
scale /i 2 = 4pf, seems to give a slightly better agreement with data than the scale /i 2 = Q 2 + p\ 
for the SaS photon structure function whereas for the DG photon structure function the scale 
/i 2 = Q 2 + p\ leads to a good description of data and fi 2 = Ap?r results in predictions which 
overshoot the data (not shown). This simply illustrates the uncertainty in the scale and in the 
parameterization of the photon structure function. The hard process q 7 g p — > qg contributes 
most to the forward jet cross section. 
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Figure 3: The ratio of cross sections for the production of two and one forward jets versus x, 
as predicted by the RAPGAP Monte Carlo for the sum of direct and resolved processes (solid 
line). The cuts applied are specified in the text. The prediction from the analytical calculation 
of [ |23| is shown with the dashed line. 

A small fraction of the DIS events, fulfilling the selection criteria for forward jets, actually 
contains two identified jets. The ratio of the cross sections for events with two and one forward 
jet(s) would be a measure of the BFKL vertex function which controls the emission of the second 
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jet. Analytic calculations (in LLA) [ ^] have been performed in the same kinematic region and 
with the same jet selection as defined above for the one-jet sample. The predicted ratio varies 
from 3% to 6%, for a cone radius of R = 1, when going from x = 0.5 • 10~ 3 to x = 3 • 10~ 3 . 
Comparisons to experimental data are not yet available. The prediction of the RAPGAP Monte 
Carlo for direct plus resolved processes is shown in Fig. 0. It is about a factor of 3 lower than 
the prediction from the BFKL calculations. A large part of this discrepancy could be due to 
hadronization effects which would reduce the prediction of the parton level BFKL calculation. 
Given the large uncertainties both in the resolved photon and BFKL calculations, it is not 
obvious that the two approaches give results in the same order of magnitude. 

4 Other tests of low x dynamics 

In this section we discuss measurements like the (2 + 1) jet rate in DIS, transverse energy 
flow and the pt spectra of charged particles, which are not well described by DIR models and 
have been subject to speculations on their possible interpretation in terms of small x BFKL 
dynamics. 

4.1 Di-jet rates in DIS 

New HI data on fractional di-jet rates were recently made public [ [|. Events with two jets 
were selected in the kinematic region 5 ~ Q 2 ~ 100 GeV 2 and 10" 4 ~ x ~ 10 2 . The jets 
were defined in the hadronic center-of-mass system using a cone algorithm of radius R — 1 and 
requiring a transverse momentum pj? > 5 GeV. The results make evident that standard DIR 
LO Monte Carlo models such as RAPGAP 2.06 in the direct mode and LEPTO 6.5 badly fail 
to describe the data. Especially in the Q 2 region below 30 GeV 2 the discrepancies are large 
but in this region p\ > Q 2 and therefore resolved photon interactions are expected to play a 
significant role. The CDM, on the other hand, gives very good agreement with data. 

With the same cuts as in the HI analysis, RAPGAP has been used, with and without 
resolved photon processes included, to predict the di-jet rates. 

The data from HI are shown as a function of Q 2 and x in Figs. |] a) and b), respectively, to- 
gether with the predictions from RAPGAP. The RAPGAP results on direct and direct +resolved 
processes are shown separately for the two scales mentioned above and for the two parameter- 
izations of the photon structure. The data are well above the direct contribution, but adding 
the direct and resolved contributions gives a remarkably good description of data all the way 
down to the phase space region where the resolved processes are dominating, especially with 
the scale fj, 2 = kp\ {ji 2 = Q 2 +p%) and the SaS (DG) parameterization of the photon structure. 
Also in this data sample the major contribution is coming from the q^g v — > qg subprocess. 
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Figure 4: The di - jet ratio -R2+1 as a function of Q 2 (a. and c.) and x {b. and d.). The points 
are preliminary HI data [|(J. In a. and b. are shown the RAPGAP predictions for the sum of 
direct and resolved processes using the scales jj 2 = Q 2 + p\ (solid line) and fi 2 = 4p|, (dotted 
line), as well as the predictions for direct processes with the scales /x 2 = Q 2 +p% (dashed line) 
and [X 2 = 4py (dashed-dotted line). The SaS 2 DIS parameterization was used to describe the 
parton density in the virtual photon. In c. and d. the sum of direct and resolved processes 
(solid line) and the direct (dashed line) are presented using the DG description of the parton 
density in the photon with the scale fi 2 = Q 2 + p 2 -. The parton density of the proton is given 
in all cases by the CTEQ4M parameterization. 
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4.2 Inclusive jet production 



The HI experiment has measured the single inclusive jet cross section in a range of photon 
virtualities from Q 2 = to 50 GeV 2 using the cluster algorithm in the hadronic CMS [ 
17). The jets from DIS processes were required to have transverse momenta above 4 GeV/c. 
Fig. [5] shows the inclusive jet cross section as a function of Q 2 in bins of E? together with the 
predictions of the DIR model and DIR+RES model. In the region where Q 2 ~ E\ the DIR 
model is expected to be valid since the photon can not be resolved. It is found that the DIR 
model approaches the data as Q 2 increases whereas in the low Q 2 region the deviations become 
significant. The DIR+RES model gives reasonable agreement with data in the full kinematic 
range investigated. 



4.3 Transverse energy flow 

The non ordered fcy parton emissions in low x BFKL dynamics produce more transverse energy 
than the standard DGLAP evolution scheme in the central rapidity region in the hadronic 
CMS, which corresponds to the forward region in the laboratory system (6 pa 10°). Results 
from HI [ H on the transverse energy flow, presented in the hadronic CMS as a function of the 
pseudo-rapidity, rj*, are shown in Fig. ^| for the kinematic region 3 GeV 2 < Q 2 < 70 GeV 2 and 
8- 10" 5 < x < 7- 10" 3 . 

The largest difference between LLA BFKL and standard DGLAP calculations are expected 
at small x and central rapidities rj*. Discrepancies between data and the DIR model of refs. [ 
are also observed in this kinematic region. Only after including a proper treatment of the 



parton dynamics of the proton remnant in the sea quark scattering processes [ and using 
a rather low cutoff on the gluon emission for the parton cascade, a reasonable agreement is 
achieved with the exception of the lowest x and Q 2 region (see Fig. However the description 
is significantly improved by including the contribution of resolved virtual photons giving an 
excellent description of the transverse energy flow over the full range of x and Q 2 as illustrated 
in Fig. H Different scales fi 2 or the parton distribution functions give changes of the order of 
5%. 



4.4 pt spectra of charged particles 



Studies based on QCD models have demonstrated that the high pt tail of charged particle 
transverse momentum spectra is sensitive to small x dynamics of parton radiation and that the 
influence from hadronization is small [ |25|, [26|]. The LLA BFKL dynamics result in harder 



spectra than obtained in the DGLAP scenario because of the non-ordered hr emissions. A LLA 
BFKL based calculation [ I27J including an estimation of the fragmentation effects was able to 



reproduce the high pr tail in the lowest x bins. 

Fig. [7| shows the pr distributions of charged particles as measured by the HI collaboration [ 
HJ for DIS events with 3 GeV 2 < Q 2 < 70 GeV 2 in the rapidity range 0.5 < rj < 1.5. The 
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Figure 5: The inclusive jet cross section as a function of Q 2 for different regions of the jet E?- 
The points are published HI data [|7|], the curves represent the RAPGAP predictions for the 
sum of direct and resolved processes (solid line) and for direct processes only (dashed line). 



HI analysis [ [Til] showed that both, direct models and CDM provide good agreement with the 
data at large x, whereas with decreasing x the direct models start falling below the data. In 
Fig. [7| the dashed line shows the prediction of the DIR model and the solid line the prediction 
of the DIR+RES model as calculated with the RAPGAP Monte Carlo program. The solid line 
in Fig. |7| demonstrates a good description of the data over the full range in x and Q 2 . 
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Figure 6: The transverse energy flow as a function of 77* in the hadronic CMS for different bins 
in Q 2 and x. The data points are preliminary HI results [ |9|]. Negative 77* corresponds to the 
proton fragmentation region. The curves give the RAPGAP predictions for the sum of direct 
and resolved processes (solid line) and for the direct contribution only (dashed line). 
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Figure 7: Transverse momentum distributions of charged particles in the rapidity range 0.5 < 
77* < 1.5 in the hadronic CMS in bins of Q 2 and x as obtained by HI [ (L(|. The RAPGAP 
prediction for direct processes (dashed line) and for the sum of direct and resolved photon 
contributions (solid line) are shown for comparison. 
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5 Discussion 



Recent experimental data on forward jet production show deviations from traditional LO Monte 
Carlo models assuming directly interacting point-like photons. This has given rise to various 
speculations. Similar deviations have been observed in other measurements like the fractional 
di-jet rates, inclusive jet production, transverse energy flow and in the transverse momen- 
tum spectra of charged particles. It is tempting to assume that the observed effects could be 
explained by BFKL dynamics as firstly the data on forward jet production, which has been 
proposed as a possible probe of low x dynamics, can not be reproduced by the DGLAP evo- 
lution and secondly that the CDM MC, which produces radiation without ordering in kx as 
expected for BFKL, reproduces all data fairly well. 

In the present study it has been shown that the addition of resolved photon processes to the 
direct interactions in DIS leads to good agreement with the data, without invoking any BFKL 
parton dynamics. We have observed that the dominant contributions to the resolved photon 
processes come from order a 2 s diagrams with the hard subprocess q^g p — ► qg (see Fig. [I]). Since 
the partons which form the photon remnant per definition have smaller px than the partons 
involved in the hard scattering, a situation with non ordering is created. 

In the LO DIR model the ladder of gluon emissions is governed by DGLAP dynamics 
giving a strong ordering of kt for emissions between the photon and the proton vertex. The 
models describing resolved photon processes and BFKL dynamics are similar in the sense 
that both lead to a breaking of this ordering in k t . The BFKL picture, however, allows for 
complete dis-ordering in k t , while in the resolved photon case the DGLAP ladder is split 
into two shorter ladders, one from the hard subsystem to the proton vertex, and one to the 
photon vertex, each of them ordered in k t (see Fig. p. Only if the ladders are long enough to 
produce additional hard radiation it might be possible to separate resolved photon processes 
from processes governed by BFKL dynamics. Thus the resolved photon approach may be a 
"sufficiently good" approximation to an exact BFKL calculation and the two approaches may 
prove indistinguishable within the range of x accessible at HERA. 

It should be emphasized again that a NLO calculation assuming point-like virtual photons 
contains a significant part of what is attributed to the resolved structure of the virtual photon 
in the RES model [ O . 
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